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We study a consistent infrared and ultraviolet regularization scheme for the cosmological pertur- 
bations. The infrared divergences are cured by assuming that the Universe undergoes a transition 
between a non-singular pre-inflationary, radiation-dominated phase and a slow-roll inflationary evo- 
lution. The ultraviolet divergences are eliminated via adiabatic subtraction. A consistent regulariza- 
tion of the field fluctuations through this transition is obtained by performing a mode matching for 
both the gauge invariant Mukhanov variable and its adiabatic expansion. We show that these quan- 
tities do not generate ultraviolet divergences other than the standard ones, when evolving through 
the matching time. We also show how the de Witt-Schwinger expansion, which can be used to 
construct the counter-terms regularizing the ultraviolet divergences, ceases to be valid well before 
horizon exit of the scales of interest. Thus, such counter-terms should not be used beyond the time 
of the horizon exit and it is unlikely that the observed power spectrum is modified by adiabatic 
subtraction, as claimed in some literature. On the contrary, the infrared regularization might have 
an impact on the observed spectrum, and we briefly discuss this possibility. 



I. INTRODUCTION 

According to the inflationary paradigm, the large scale 
structure observed today in the sky originated from tiny 
quantum fluctuation in the early Universe. At that time, 
the Universe experienced a period of quasi-exponential 
expansion that amplified these fluctuations, which be- 
came classical after exiting the horizon and generated the 
gravitational instabilities responsible for the formation of 
the large-scale structures. One of the most important 
predictions of the inflationary scenario is that the spec- 
trum of these fluctuations is nearly scale- invariant pQ. 

The quantum origin of the perturbations of the inflaton 
field and of the metric tensor necessarily raises a concern 
about renormalization. In fact, it is well known that the 
correlation functions of quantum fields on a curved back- 
ground suffer from divergences that cannot be cured as in 
flat space [2] . In the case of a time-dependent Friedmann- 
Lemaitre-Robertson- Walker (FLRW) background, quan- 
tum correlations typically develop logarithmic diver- 
gences in the ultraviolet (UV), together with the famil- 
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iar quadratic ones. In addition, also infrared (IR) di- 
vergences exist. These problems become quite relevant, 
as the observed power spectrum is directly connected to 
the two-point function of the gauge invariant Mukhanov 
variable [3] in the coincidence limit. 

Concerning the UV divergences, the infinities can be 
cancelled by subtracting counter-terms constructed ac- 
cording to the adiabatic expansion in the momentum 
space [4], which is equivalent to the de Witt-Schwinger 
technique in coordinate space [5]. According to [SHE], 
the adiabatic subtraction leaves an imprint on the renor- 
malized power spectrum because the adiabatic counter- 
terms, when evaluated a few Hubble time after the hori- 
zon exit, are relatively large for the scales of interest. On 
the other hand, in a recent paper [Hj , we argued that this 
is not correct. Among other problems, the main one is 
that the adiabatic subtraction procedure seems to be ill- 
defined when the scales of interest are stretched towards 
horizon exit (see also [TU] for a different criticism on the 
main idea proposed in [5] ). At the end of this paper, 
we offer a further explanation on why the adiabatic sub- 
traction should not be considered valid around horizon 
crossing. 

Regarding the cure for IR divergences one possibility 
consists in assuming an initial vacuum state which dif- 
fers from the usual Bunch-Davies vacuum. Note that this 
approach is different from what was done in the context 
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of the trans-Planckian problem, see e. g. [TT], where the 
interest was about the UV behavior only. Physically, this 
is equivalent to assume that the Universe emerges from a 
pre-inflationary phase dominated, for example, by mat- 
ter or radiation, see e. g. [T5J [T3] [21] ■ A preinflationary 
matter-dominated phase, and its effect on the CMB, was 
also considered in [15j . Then, in line with an old theorem 
formulated by Ford and Parker (16], no IR divergences 
can develop during the subsequent expansion. In this 
paper, we assume that fluctuation modes evolve across 
a sharp transition from a radiation-dominated Universe 
to an inflationary slow- roll phase [57]. The only require- 
ment that we impose is the continuity of the scale fac- 
tor and of its first derivative as in [12] . Mode matching 
has been already utilized to study observational signa- 
tures of pre-inflationary phases characterized by lower- 
dimensional effective gravity or modified dispersion rela- 
tion in |18j . Also, the spectrum of gravitational waves 
generated by a series of radiation and matter dominated 
phases has been studied in [19] . 

The main question that we want to address here is 
whether the UV and IR renormalization schemes out- 
lined above can be performed together. Specifically, first 
we show that the new terms that arise from a non-trivial 
vacuum, and that regularize the IR divergences, do not 
generate, after the match, new UV divergences with re- 
spect to the standard ones which characterize the infla- 
tionary phase. Then, we want to verify if also the adia- 
batic counter-terms, that must be present already before 
the match, evolve consistently through the match. 

The plan of the paper is the following. In Sec. [TT] we 
introduce the formalism. In Sec. Mil we outline the mode 
matching technique applied to the scalar perturbations 
in terms of the Mukhanov variable. In Sec. |IV| we study 
the propagation of the IR-regulating terms through the 
phase transition at the matching time. In Sec. [V] we ex- 
tend this analysis to the adiabatic counter-terms. In Sec. 
|VI| we estimate the potential impact of these new terms in 
the prediction of the inflationary spectra. As mentioned 
above, in Sec. VII we show why adiabatic subtraction 



should not be taken too seriously when the scales of in- 
terest to us cross the horizon. We finally conclude in Sec. 



VIII with a summary of our results. 



the background equations of motion for 0(t) and for the 
scale factor a(t) read 
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2M p 2 , 
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where M 2 Y = l/(8nG) is the reduced Planck mass, and 
the energy density p and the pressure p arc related by the 
equation of state p = up. The dot denotes a derivative 
with respect to the cosmic time t. 

The perturbations of the background metric can be 
written in the form 



ds 2 



-(1 + 2^)dt 2 + a 2 [(1 - 2$)5ij + hij] dx l dx j , (6) 



and for single-scalar field inflationary models the Bardeen 
potentials \& and <& coincide. Analogously, also the infla- 
ton field is perturbed according to <f> — > <j> + 8(f). One can 
conveniently describe the scalar perturbations by means 
of the so-called Mukhanov variable Q, defined as [3] 



= 4- 



2H- 1 i> + * 
3(l + u;) 



(7) 
(8) 



For the second expression, which is also well defined in a 
Universe dominated by a fluid, we have used the Einstein 
constraint equation, see e.g. [50]. Upon quantization, one 
promotes the variable Q to the operator 



Q(t,x) 



1 



(2tt) 3 /- 



Qk(t)e 



ik-x 



Q* k (t)e 



ik-x it 



(9) 

where bk is a time-independent Heisenberg operator that 
satisfies the usual commutation relations 



[b k ,k'] = [i>l bl\ = 0, [k,i>U =<5 (3) (k-k') , (10) 
provided the modes Qk satisfy the Wronskian condition 



II. FORMALISM 



,3 • 



(11) 



Let us consider a spatially flat Universe, whose dynamics 
is driven by a classical minimally coupled scalar field, and 
is described by the action 



S = I d i xy/^g 



For a spatially flat FLRW spacetime, with metric 
ds 2 = -dt 2 + a 2 (t)6 ij dx i dx j , 



(1) 



(2) 



The equation of motion is then given by 



Qk + iHQk 









V 4>4> 




-?)] 









= o, 

(12) 

where denotes the second derivative of V with re- 
spect to 4>. We now introduce the slow roll parameters: 



2M*H* 



H 

IP 



V 



(13) 
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Note that our definition of e is more general than the 
usual one, e = (M 2 l /2)(V ( f,/V) 2 , as it is well-defined for 
arbitrary backgrounds, not only the ones dominated by 
a scalar field. With these, Eq. ( 12 ) can be rewritten as 



k 2 

Q k + 3HQ k + ^Q k + 
a z 



3r? - 6e + 2e 2 



rje 



H 



Qk=0- (14) 



We note that the derivative e is related to e and rj by 
e = -2e(?7 - 2e)H . (15) 



This equation can be used to replace e in Eq. ( 14 ) but it 



can also be used as a definition of rj for an arbitrary mat- 
ter content of the Universe, not necessarily a scalar field. 



Eq. ( 14 1 holds only if the degree of freedom Q comes from 



a scalar field. If it corresponds to fluctuations in a fluid 
with equation of state of the form p = ojp, the k 2 term 
must be multiplied by the factor c 2 = p/p. In all other 
aspects, the linear perturbation equation remains identi- 
cal, see [2U] for more details. Finally, Q k is associated to 
the scalar power spectrum, defined by 



P c {k) 



k 3 ^ 

2TT 2 



±-K 2 M^6 



\Q k \ 



(16) 



where we used Eq. ( 13 ) for the second equality. 



We now consider two different cases of interest for our 
investigation. The first one is a slow-rolling inflationary 
evol utio n where e <C 1 and where, as one can see from 
Eq. ( 15 ), e = 0(e 2 ) can be neglected to the leading order. 



The second one is the case when e is exactly constant 
so that rj = 2e, again because of Eq. (15). The latter 
describes a power law expansion, as for the case of a 
radiation or matter-dominated Universe. In both cases 



the general solution to Eq. ( 14 ) can be written in terms 
of Hankel functions, namely 



Ik = a 



[E(k)u(z) + F{k)u*(z)} 



where 



u(z) 



-H^(z) 



and where we have introduced the new variable 
z(t,k) 



(1 - e)aH 



(17) 



(18) 



(19) 



The index of the Hankel function is given by (we neglect 
contributions which are subleading for the slow-roll in- 
flationary case or exactly zero for the case e = 0) 



2(1 -e) 



3ry + 6e . 



(20) 



The Wronskian condition (11) implies that |.E(fc)| z 
\F(k)\ 2 = l. 



When e is constant, it is related to the equation of state 
via e = |(l+w). In Fig. IT] we plot e and fi as functions of 
lu. Note the divergence at w — —1/3, which corresponds 
to a curvature dominated universe. In this case e = 1 and 
z is not well defined. Note also that, since e = 1 — ad/ a 2 , 
we have that e < 1 for an accelerating Universe and e > 1 
for an decelerating one. Typical cases are 

• Radiation-dominated: u) = 1/3 and e = 2, fj, = 1/2. 

• Matter-dominated: u> = and e = 3/2, /x = 3/2. 

• A-dominated: in the slow-roll approximation we 
have typically, at the leading order, /i = 3/2 — 77 + 3e 
with < e < 1 and < \rj\ < 1. Thus the Uni- 
verse accelerates. For the exact de Sitter solution, 
one has ui = — 1, e = rj = and fx = 3/2. 

The first two cases can be obtained either by a fluid 
with the corresponding equation of state or by a scalar 
field with exponential potential, which yields a power-law 
scale factor (see, for example, [2"T]). 
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FIG. 1: We plot e (red dashed line) and \fi\ (blue solid line) 
as functions of w for power law expansion. 

We remark that (only in four space-time dimensions) 
fx = 3/2 represents a "degenerate" solution, as it de- 
scribes both an accelerating, de Sitter Universe, and a 
decelerating, matter-dominated Universe. In this case, 
one must specify also u or e. Note also that, when the 
expansion of the Universe decelerates one obtains z < 0. 
Therefore, the solution ( 18 ) should be modified to 



(-*) 



(21) 



as one usually takes the negative axis as the branch cut 
for Hankel functions and one has to satisfy the Wronskian 
condition. 

During inflation, the wavelength of a mode k is grow- 
ing, A = (2w/k)a(t), while the Hubble parameter remains 
nearly constant, hence z is decreasing. A typical mode 
is inside the horizon at early times, z » 1 and crosses 
the horizon as inflation unfolds. For such a mode we 
can argue that it is initially in the Minkowski vacuum, 
u(z) oc exp(ikrj) / V2~k = exp(—iz)/\2k. This leads to 
the inflationary initial condition given by Qk{z) — 
= 1 and F„ = 0. 
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However, during a decelerated expansion, \z\ is grow- 
ing and modes which are inside the horizon when inflation 
begins, have been outside during the decelerated expan- 
sion before inflation, and in general we have no way to 
determine their initial condition. For large scales, time 
dependence cannot be neglected and for general time de- 
pendent spacetimes we cannot formulate vacuum initial 
conditions. 

Fortunately there is one exception to this rule and this 
is the radiation dominated Universe: in this case, we can 
redefine our perturbation variable to 



v = aQ 



(22) 



such that Eq. (14 1, in terms of conformal time t\ defined 
by dr\ = dt / a(t)^reduces to the simple Minkowski wave 
equation 



k 2 v = 



(23) 



Here a prime denotes the derivative w.r.t. conformal 
time 77. (Again, if we replace the scalar field lead- 
ing to an radiation-dominated expansion law, a oc i 1 / 2 , 
with a radiation fluid, we must substitute the term k 2 
by k 2 /3.) As we know how to quantize Eq. (23 1, we 



can set up Minkowski vacuum initial conditions for all 
modes. Here the expansion of the universe has disap- 
peared and is taken into account simply by the normaliza- 
tion, v — aQ and dt = adrj. The vacuum initial condition 
for v corresponds exactly to Eq. (21 1 for n = 1/2, and 



v = aQ = u(z) = ^/^H^li-z) = iexp(-ikr])/V2k. 
Apart from the de Sitter case, where time dependence 
can be fully removed, this is the only FLRW background 
which allows for well defined quantum initial conditions 
for all modes. The reason for that is that for a radia- 
tion dominated universe, the Ricci scalar, R = H + 2H 2 
vanishes. Therefore, any scalar field can be considered 
as conformally coupled and as such is not affected by the 
expansion of the Universe. 

This observation is crucial for the model which we de- 
velop in the next section where we match a previous ra- 
diation era to a subsequent inflationary era. 



III. MODE MATCHING AS A CURE FOR IR 
DIVERGENCES 

We now assume that the solution to the mode equation 
is characterized by two cosmological phases determined 
by two different values of the index [i, which we denote 
by jjL and a. The global evolution of Q then is 



(24) 



Q a = a-\t) [E a (k)u a (z) + F a (k)u* a (z)} , U < t. (25) 



As we mentioned above, Eq. ([24| is a sensible vacuum ini- 
tial condition only for a radiation-dominated or inflation- 
ary phase. We shall be interested in the first case later 
on, but we keep our discussion general in the beginning. 



We impose that a^U) = a a (ti) and d M (ij) = a a (ti). 
This can be realized in the following way. Consider the 
equation e = —H/H 2 = const. The general solution for 
the Hubble parameter is 



H = 



Ho 



l + eH {t-t o y 
which implies that 

o(t) = a [l + eH (t-t )] 1/e . 



(26) 



(27) 



In these expressions, Hq and aQ denote the arbitrary val- 
ues of H (t) and a(t) at t = t . Now, suppose that 



a(t) 



fll [l +£!#!(< -t!)] 1/£l t<t t 

a 2 [l + e 2 H 2 (t-t 2 )} 1/e2 t>U 



(28) 



If we consider the case for which Hi = H 2 the match can 
be easily imposed setting ti = t 2 = ti and a\ = a 2 . On 
the contrary, if Hi 7^ H 2 , by imposing continuity of H (t) 
across ti we find 



U = 



H 2 — Hi 



HxH 2 {e 2 -ex) ' 
together with the consistency condition 
eiti = e 2 t 2 = c , 



(29) 



(30) 



where c is an arbitrary constant. The continuity of a(t) 
across ti further implies that 



«2 

a 1 



H 2 e 2 — Hi€i 
£2 - £l 



cHiH 2 



-| l/ei-l/e 2 Jf 1 /^ 



H 



(31) 



We are interested in a transition between a radiation- 
dominated Universe (ei = 2) and a slow-roll inflation- 
ary one (0 < e 2 <C 1). In this case, the scale fac- 
tor has the form a(t) ~ t 1 ! 2 before the match and 

a(t) ~ aoe^^'o'+Trt*-' ' 2 after the match. So con- 
sidering an inflationary phase that starts at the time of 
the match the scale factor takes the following form 



a(t) 



{t/kf/ 2 t<t 



(32) 



t > k 



with Hi = 277. The inflationary expression for the scale 
factor is exact for a quadratic potential, V — m 2 (j) 2 /2, 
and is a first order expansion in the slow-roll parame- 
ters for other choices of the inflationary potential. For 
the numerical results shown below we always choose a 
quadratic potential. 

We now determine the coefficients E n and F n in 



Eq. ( 25 ) , by solving the linear system 

Q li(j^i) Qa{pi) 
Q ^li^i) — Qaiti) 



(33) 
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with the help of the result 

m * - 



k 



(34) 



which holds whenever e is nearly constant, and using the 
Wronskian identity (7 = a, /x) 



we find 



E a = ia(ti) (u* u M - u* a u^) t=ti , 
F Q = ia(t 4 ) (u Q u p - w Q u M ) t=t . , 



(35) 



(36) 



and it is easy to check that | E a | 2 — | F a \ 2 = 1 . 

In the rest of the paper, we assume that the state of 
the Universe preceding the inflation is dominated by a 
scalar field with an exponential potential leading to a 
radiation-like expansion. To realize this, it is sufficient 
to take ED 



V 



Mpi In [ - 

K 

4t? 



2m) - 4>i 



M, 



pi 



(37) 



where <pi is an arbitrary constant. In fact, one can check 
that e = 2, 77 = 4 thus, from Eq. ( 14 1 one obtains 



Q = 



2k 



(38) 



Another advantage of this point of view is that the 
Mukhanov formalism is well-defined for all times, and 



there is no need to adjust Eq. ( 14 ) with the speed of 
sound. 

The infrared divergences become apparent when one 
calculates the two-point function G at coincident points 



Goc J dkk 2 \Q c 



In the small z limit, one finds that 

u a (z) ■ 




u* a (z) 



(39) 

(40) 
(41) 



thus |Q Q | 2 cx \u a \ 2 \E a -F a \ 2 ex k- 2a \E a ~F a \ 2 . By using 
Eq. ( |40[ ), its time derivative and similar expressions for 
u M one finds that, for z(U) < I, \E a ~ F a \ 2 ~ k 2 ^'^ 
and 



G cx / dkk 2{1 -^ . 



(42) 



This result first shows that the low-fc behavior of the inte- 
gral is independent of a, namely of whether the Universe 
accelerates or decelerates after the transition at t = ti- 
Then, it is clear that the integral converges provided 
/i < 3/2, which excludes configurations where the Uni- 
verse accelerates for t <U, such as de Sitter or slow-roll. 
On the contrary, if the Universe begins in a radiation- 
dominated phase /i = 1/2 the IR convergence is guaran- 
teed. These calculations confirm the statement that IR 
divergences cannot, in general, develop during a smooth 
evolution of the Universe |16j . Similar results were found 
in [T3] , where the evolution of a test scalar field was stud- 
ied during a smooth transition from a decelerating and 
expanding Universe to an inflationary one. 



IV. UV IMPACT OF THE MODE MATCHING 

The above results show that mode matching can solve the 
IR problem if the very first phase of the Universe is of 
non-inflationary type. The infrared end of the spectrum 
is sensible only to this phase, and cannot be changed by 
the future evolution (provided e is nearly constant). On 
the ultraviolet side, we know that adiabatic subtraction 
of appropriate terms can cancel the UV divergences. The 
main problem, which has not been addressed yet, is the 
evolution of these counter-terms through the transition 
at t = ti. We will address this issue in the next section, 
but first let us see if the divergent structure after the 
match is influenced by the presence of the match itself. 

After the transition, the wave function is no longer in 
a Bunch-Davies state, as it is represented by the function 



Q a = [E a (k)u a (z) + F a {k)ul{z)} 



(43) 



where E a and F a are given by Eq. (36). We wish to 
evaluate \Q a \ 2 in the ultraviolet regime. For large z we 
have [22] 



[Sf,, + iRn] , 



(44) 



where 



Rl 



S, = 



(4 M 2 -l)(4^-9) 
128z 2 



4 M 2 - 1 

> 1 

Z- 7T I - + - 

2 4 

r(2) 



(45) 

(46) 



and similar expressions for H£ J . With these expansions 
and Eq. ( 36 ) , calling z M the value of z for t < and z a 
the value for t > ti, we find that 



\Qc 




1 



2fca 2 ) 8zl 
4a 2 - 1 4^i 2 - 1] cos(2z Q - 2z a ) 



(47) 
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where the barred quantities are those evaluated at t = U. 
The third term in square brackets only appears because 
of the match at t — ti, and consistently vanishes when 
fi = a. This term is not divergent in the UV for t > U. 
Before the match, the ultraviolet structure of the modes 



is 



1 



(V - 1) 



(48) 



and shows the usual logarithmic and quadratic diver- 
gences (except for the radiation-dominated case consid- 
ered, where = 1/2 and the logarithmic one is not 
present). From Eq. (47) we see that the divergent struc- 



ture after the match is not altered by the presence of the 
match at t = ti. As a check of the method, one easily ver- 
ifies that at t = ti the ultraviolet limit of |Q Q | 2 coincides 
with the one of IQ^I 2 , as required by the first condition 
in Eqs. (pi 



V. 



ADIABATIC EXPANSION THROUGH THE 
MATCH 



The UV regularization of the two-point function can be 
achieved by subtracting appropriate counter-terms from 
the divergent integral. These quantities can be con- 
structed by solving, through an adiabatic expansion at 
the appropriate order, the mode equation [51 0] . In order 
to be consistent with the mode matching, we wish com- 
pute the counter-terms and study their evolution through 
the phase transition at t = ti. Let us write the adiabatic 
solution, before and after the match, as 



We stress that also the term 
of second order as well l28l. 



should be considered 



By solving again the system ( 33 1 for the adiabatic case, 
we find the coefficients 



pad 



(ilka + Qkfi) 



2(fi fcQ f>^)i/2 



(53) 
(54) 



where only terms up to the second adiabatic order have 
been retained. With these, we find 



ad|2 



1 



4a2f2 fcQ 
x cos I 2 



dt 1 



^ka ^kfi 
*^kfj, ^ka 



-Qkr 



(55) 



Expanding up to the second adiabatic order this expres- 
sion can be written as 



ad 1 2 



1 



2alk 
cos I 2 



2 k 2 



Ml) 



dt' 



t, a a (t') 



f \ ^ k ot 



(56) 



The first part of this expansion cancel the divergent terms 
in Eq. (47) in the ultraviolet limit. So, using such expres- 



sion we can have a consistent regularization of the field 
fluctuations through the time of the matching and be- 
yond. 



l 



a, 



,^2~W } 



: CXp I 



kfi. 



dt. 



t 



-w k . 



Q a d (t) = 4r exp ! / - w *< 

a a y/2W ka \ J t . a a 

pad / ft j. \ 

+ /^ exp H / ~ Wk A > ( 49 ) 

a a ^/2W ka \ J u a a J 

where the superscript "ad" stands for adiabatic, and 
where W ks must satisfy the equation 



w 2 kl 



«* 7 



1 



w kl 



(50) 



where 7 = [if a. The generalized frequency il kl is defined 
as 



n 2 kl = k 2 + a 2 M 2 



with 



M 2 = V 



R-6 



2a" 



,'2 



2a 



112 



(51) 
(52) 



where the underscript 7 is neglected in the r.h.s.. For our 
purposes, it is sufficient to keep the terms up to the sec- 
ond adiabatic order (namely, with two time derivatives). 



VI. OBSERVATIONAL SIGNATURES 

We want now to evaluate possible signatures of the match 
considered on the spectrum of the scalar perturbations. 
Every mode exiting the horizon at the time t cx satisfies 
the relation a(t ex )H(t cx ) — k. In typical inflationary 
models, modes that exit about 60 e-folds before the end 
inflation correspond to scales that are observable today. 
We wish to compute the spectrum associated with these 
modes and verify whether the pre-match phase can have 
left some signature. Thus, we calculate the coefficients 
E a and F a at the time t = U according to Eqs. (36) and 
|Q Q | 2 according to Eq. (43). Then we can compute the 



scalar power spectrum defined by Eq.([16j). 

Let us consider two limiting cases, the one where U — 
t cx and, as consequences, a(ti)H(ti) = a(t cx )H(t cx ) = 
k and the one where ti -C t cx with a(ti)H(ti) <C 
a(t cx )H(t cx ) 

using Eqs. ( 44|45 ), one finds 



k. In this second case z^,z a 3> 1 and, 



\E n 



l + O 



1 



\F„ 



O 



1 



E a F* = O 



(57) 



(58) 
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At the leading order, we have 



\Q a \ 



1 



[(\E a \ 2 + \F a \ 2 ) 



EaKK 



(59) 



so we recover the standard slow-roll behavior, and we 
have no observational consequences of the match at the 
leading order. 

On the contrary, in the case U = t cx , the contribu- 
tion of E a and F a give sensible corrections to the power 
spectrum. To see this, let us first evaluate the spectrum 
exactly at the horizon exit. In the standard case we have 



(60) 



with z a (k — aH) ~ (1 + e). By expanding with respect 
to the slow-roll parameters, the leading term is given by 
the Hankel function with a = 3/2, so 



and 



P ( st (k = aH) 



2 

TTZ n 



1 



1 



(61) 



(62) 



On the other hand, if one considers the matching condi- 
tion Q a (ti) — Qn(ti) together with 



1 I * u m ( — 
AaH " \aH 



(63) 



with (i = 1/2 (radiation domination), the modified spec- 
trum at the horizon exit becomes 



P c (k = aH) 



1 



H 



2 .\ /,-;.< V 2tt 



(64) 



so it is reduced by about 50%. 

Let us now instead consider the case in which the spec- 
trum is evaluated several e-folds (already 5 e-folds are 
sufficient) after the horizon exit, namely when k <C aH. 
This means that the argument of the Hankel function 
is relatively small and we can make the approximation 

l-Ho^C^a)! 2 ~ 2/(wz^ l ). Then, the spectrum has the well- 
known value 



P c st {k < aH) 



1 

2M 2 



(65) 



As before, the introduction of the matching conditions 
together with the requirement that the fluctuations exit 
the horizon at the time of the match (k = ciiHi) changes 
our result. In the limit k <C aH we have that it* ~ — u a 
and, as in Sec. Ill, one obtains 



\Qc 



I -En 



(66) 



The first part on the r.h.s. gives the standard contri- 
bution to the spectrum, while \E a — F a \ 2 calculated on 
k = aiHi is a numerical coefficient, which is independent 
of the initial condition and nearly equal to 0.38. So the 
spectrum, in the presence of the match, is reduced by 
about the 62%. 

The general behavior of the modified spectrum for 
these two cases is shown in Figs. [2] and [3] Here, we plot 
the ratio of the modified power spectrum with respect 
to the standard one for the case when the inflationary 
potential is V = ^m 2 (f) 2 , with the match described in 
section III and an initial condition fixed by a value of 
H(ti), which guarantees near 60 e-folds of inflation. 




FIG. 2: P c (k = aH)/P ( at (k = ah) is shown for a mV po- 
tential versus the number of e-folds iV between the beginning 
of the inflation and the time where the fluctuation considered 
crosses the horizon. The spectra are calculated at the time of 
horizon exit. 



To conclude this section, let us consider fluctuations 
with wave numbers k < a(ti)H(tj), namely fluctuations 
that were outside the horizon at the beginning of infla- 
tion, and still are today. Although their spectrum is not 
directly accessible to observations, these modes can have 
an effect on non-Gaussianities produced at second or- 
der, which will be the subject of a future work. If wc 
consider the limit k <C a{ti)H{ti) < a(t)H(t) we have 

the same condition that leads to Eq. (66) but now, as 

k 2( a - 



shown in Sec. Ill, \E a — F a \ 2 ~ k 2 ^ a ^~and the spec- 
trum changes from nearly scale-invariant to very blue, 
P c (k) ~ k 2 . In Fig. g we show the general behavior of 
this spectrum, calculated several e-folds after the exit of 
the observable fluctuations, going from k -C a(ti)H(ti) to 
k ^> a(ti)H(ti). As discussed in [9], the renormalization 
of the spectrum by adiabatic subtraction at the horizon 
exit seems meaningless (see also next section). Therefore, 
we omit the evaluation at horizon exit of the adiabatic 
counter-terms. 
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FIG. 3: P c (k < aH)/P ( st (k < aH) is shown for a m 2 cj> 2 po- 
tential versus the number of e-folds N between the beginning 
of the inflation and the time where the fluctuation considered 
crosses the horizon. The spectra are calculated several e-folds 
after the horizon exit. 




FIG. 4: log (P ( (k<^aH)/P ( st (k<^aH)) versus the 
log(fc/m) is shown for a m 2 <f> 2 model with nearly 60 e-folds 
of inflation. The range of k goes from 0.1m to the value that 
exits the horizon 3 e-folds after the beginning of the inflation 
(for k = a(ti)H(ti) we have log(fc/m) ~ 1.85). The spectra 
are calculated several e-folds after the beginning of inflation. 



VII. VALIDITY RANGE OF THE ADIABATIC 
SUBTRACTION 

In this section, we would like to study the validity range 
of the adiabatic expansion. The adiabatic counter-terms 
obtained with an adiabatic expansion, as the one shown 



in Sec. [Vl can be found also with a de Witt-Schwinger 
(dWS) point-spitting. However, we shall see that the 
dWS series is no longer valid in the slow-roll approxima- 
tion as one approaches horizon exit. From the point of 
view of the adiabatic subtraction, this is equivalent to 
the loss of adiabaticity. 

To show this, let us briefly recall the dWS procedure 
(for more details, see [HE]). The starting point is the 
local expansion of the metric in Ricmann Normal Coor- 
dinates (RNC), which can be regarded as a constructive 
proof of the local flatness theorem [2J [23| ■ On a smooth 
manifold, we can always expand the metric around a 
given point P as 



(67) 



where y a are the RNC with origin at P, rjuv is the 
Minkowski metric, and the dots represent higher curva- 
ture terms. If the vector with components p a is the tan- 
gent at P to the geodesies that joins the points P and 
Q, we can define the RNC as y a = rp a , where r is an 
affine parameter along the geodesies. It is clear, from the 
construction itself, that the validity of the RNC patch is 
restricted to a region where geodesies do not intersect. 
There is however another constraint on the typical size 
L of the volume around P where RNC are valid, namely 
that the curvature is not rapidly changing. These validity 
regimes can be stated more precisely as jMj 



L <C min 



1 



\R 



(68) 



where the first bound comes from the size of volume with- 
out intersections and the second from the size of volume 
where the curvature is slowly varying. If these condi- 
tions are met, one can define a local Fourier transform 
operator and write the two-point function as [5] 



G(P,Q) 



d 4 k 



iky G{k) 



(69) 



For a minimally coupled massless scalar field, the func- 
tion G(k) can be expanded as 



afK 

6p~ 



+ 



(70) 



with the expansion valid only when \a 2 R/k 2 \ <C 1. On a 
spatially flat FLRW background, this condition reads 



a 2 H 2 
k 2 



H 

H 2 



< 1 



(71) 



In the particular case of slow-roll inflation we see that 
the validity of the RNC expansion is set by the constraint 
aH <C k, so it is not correct to calculate spectra at the 
horizon exit, aH = k, with this renormalization method. 
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We can verify this finding by looking directly at the 
bounds imposed by the inequalities (68). For a FLRW 
background with metric ^ one has 



RoiOj — (e — l)a 2 H 2 8i 



D _ „4 rr2 



(72) 



where latin indices label spatial coordinates and there is 
no summation over repeated indices. It follows that, in 
the slow- roll regime, when e ~ and |e| <C 1, one has 



R 



OiOj 



r 



OiOj 



1 



1-eH 



Ri 



Ri 



1 



2 2-eLff 



(73) 



and this implies that the dWS expansion is valid in a 
region of size L <C l/(aH), namely much smaller than the 
comoving Hubble radius that decreases during inflation. 
To enlarge the validity domain, one should consider more 
terms in the expansion. But this would spoil the rule that 
the number of adiabatic counter-terms has to be two, 
in order to just cure the quadratic and the logarithmic 
divergence. 

These considerations reinforce the findings of [S], 
namely that adiabatic or dWS subtraction are not suit- 
able to renormalizc the two-point function, and hence the 
power spectrum, at and after the horizon exit. 



VIII. CONCLUSIONS 

In this paper we have examined the simultaneous reg- 
ularization of both infrared and ultraviolet divergences 
of cosmological perturbations, treated as quantum fields 
on a curved background. The infrared regularization is 
performed by matching slow roll inflation to a radiation 
dominated, pre-inflationary phase. We do not take this 
particular phase too serious, but our procedure shows, 
with this concrete example, that a previous phase with a 
well defined initial vacuum can well regularize inflation- 
ary infrared divergences. It will be interesting to study 



in the future, whether higher order perturbations, which 
are relevant e.g. for non-Gaussianities, will depend on 
the details of this infrared regularization. Furthermore, 
since our modification happens before the onset of infla- 
tion and just removes power in the modes which never 
enter the horizon during inflation, we expect it to regular- 
ize the infrared also for inflationary models which deviate 
from slow roll as, for example, warm inflation [25 . The 
details of the matching and the possible observational 
signatures might however be somewhat modified. 

To regularize the ultraviolet divergence we employ the 
usual subtraction of adiabatic counter-terms. The first 
result is that the mode matching does not introduce new 
ultraviolet divergences. In addition, the adiabatic expan- 
sion is well defined through the match despite what one 
might fear because of the discontinuity of the Ricci scalar 
R. In fact, in our model R is not continuous through the 
transition as it contains second derivatives of the scale 
factor. Therefore, one might expect that the de Witt- 
Schwinger expansion ( |70[ ) , which is equivalent to the adi- 
abatic expansion, be no longer well defined. On the con- 
trary, we find that this is not the case. 

Finally, we argue that the adiabatic expansion is not 
valid at and after the horizon exit, by looking more care- 
fully at the construction of the counter-terms by the de 
Witt-Schwinger point-splitting method. This reinforces 
our opinion that ultraviolet regularization cannot leave 
any observable imprint. However, it seems possible that 
the mode matching in the infrared has left some obser- 
vational trace, but only if it occurs very close to horizon 
exit of the scales of interest. 
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